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Abstract 

The relativistic mean- field plus RPA calculations, based on effective Lagrangians with density- 
dependent meson-nucleon vertex functions, are employed in a microscopic analysis of the nuclear 
matter compressibility and symmetry energy. We compute the isoscalar monopole and the isovector 
dipole response of ■^''^Pb, as well as the differences between the neutron and proton radii for ■^''^Pb 
and several Sn isotopes. The comparison of the calculated excitation energies with the experimental 
data on the giant monopole resonance in ^'^^Pb, restricts the nuclear matter compression modulus of 
structure models based on the relativistic mean-field approximation to ~ 250 — 270 MeV. The 
isovector giant dipole resonance in ■^''^Pb, and the available data on differences between neutron 
and proton radii, limit the range of the nuclear matter symmetry energy at saturation (volume 
asymmetry) to 32 MeV < 04 < 36 MeV. 

PACS numbers: 21.60.Ev, 21.60.Jz, 21.65. +f, 24.30. Cz 
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I. INTRODUCTION 



Basic properties of nuclear ground states, excitation energies of giant monopole reso- 
nances, the structure of neutron stars, the dynamics of heavy-ion colhsions and of supernovae 
explosions, depend on the nuclear matter compressibility. The nuclear matter compression 
modulus Knm is defined as 



_ d-E/A 



(1) 



where E/A is the binding energy per nucleon, kf is the Fermi momentum, and kf^ is the 
equilibrium Fermi momentum. The value of -ft'nm cannot be measured directly. In principle 
it can be extracted from the experimental energies of isoscalar monopole vibrations (giant 
monopole resonances GMR) in nuclei. Semi-empirical macroscopic leptodermous expansions, 
as well as microscopic calculations, have been employed in the analysis of available data 
on isoscalar GMR. Although macroscopic expansions, analogous to the liquid drop mass 
formula, in principle provide "model independent" estimates of K^^, in reality they do not 
constrain its value to better than 50%. 

A more reliable approach to the determination of is based on microscopic calcula- 
tions of GMR excitation energies. Self-consistent mean-field calculations of nuclear ground 
state properties are performed by using effective interactions with different values of -ft^nm- 
Interactions that differ in their prediction of the nuclear matter compressibility, but other- 
wise reproduce experimental data on ground-state properties reasonably well, are then used 
to calculate GMR in the random phase approximation or the time-dependent framework. A 
fully self-consistent calculation of both ground-state properties and GMR excitation ener- 
gies restricts the range of possible values for K^rn- The correct value of K^^n should then be 
given by that interaction which reproduces the excitation energies of GMR in finite nuclei. 
It has been pointed out, however, that, since -f^nm determines bulk properties of nuclei and, 
on the other hand, the GMR excitation energies depend also on the surface compressibility, 
the determination of the nuclear matter compressibility relies more on a good measurement 
and microscopic calculations of GMR in a single heavy nucleus such as ^°^Pb, rather than 
on the systematics over the whole periodic table . It has also been emphasized that the 
determination of a static quantity, -ft'nm, from dynamical properties, i.e. from GMR energies, 
is potentially ambiguous. Various dynamical effects, as for instance the coupling of single 
particle and collective degrees of freedom, could modify, though not much, the deduced value 
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of the nuclear matter compression modulus. On the other hand, it has been shown that K„ m 
cannot be extracted from static properties, i.e. masses and charge distributions, alone jsf. 

In this work we address a different source of ambiguity, which has become apparent only 
recently. Modern non-relativistic Hartree-Fock plus random phase approximation (RPA) 
calculations, using both Skyrme and Gogny effective interactions, indicate that the value 
of -ft'nm should bc in the range 210-220 MeV Q, Q. In Ref. Q it has been shown that 
even generalized Skyrme forces, with both density- and momentum-dependent terms, can 
only reproduce the measured breathing mode energies for values of i^nm in the interval 
215 ± 15 MeV. A comparison of the most recent data on the EO strength distributions 
in ^^Zt, ^^^Sn, ^^''Sm and ^°^Pb, with microscopic calculations based on Gogny effective 
interactions by Blaizot et al. jl|, has put the value of K^^^ at 231 ±5 MeV J^. In relativistic 
mean-field (RMF) models on the other hand, results of both time-dependent and RPA 
calculations suggest that empirical GMR energies are best reproduced by an effective force 
with JsTnin ~ 250 - 270 MeV 

mm 

. Twenty percent, of course, represents a rather large 
difference. The origin of this discrepancy is at present not understood. 

It appears, however, that not all results of the relativistic models are in agreement either. 
By using the same type of non-linear RMF Lagrangians with scalar self-interactions as 
in Refs. series of recent papers j9|, |10|, |ll| Piekarewicz has performed 

fully consistent relativistic RPA calculations of nuclear compressional modes. Basically the 
results can be summarized as follows. (1) The best description of both compressional modes 
- isoscalar monopole and dipole - is obtained by using a "soft" effective interaction with 
i^'nm = 224 MeV. (2) Effective interactions with compression moduli well above -ft'nm ~ 200 
MeV are not consistent with experimental data. (3) The compression modulus determined 
from the empirical excitation energies of the GMR depends on the nuclear matter symmetry 
energy, i.e. models with a lower symmetry energy at saturation density reproduce the 
GMR in ^''^Pb by using a lower value of -ft'nm- (4) The variance between the values of 
fCnm determined from non-relativistic and relativistic mean-field plus RPA calculations of 
GMR excitation energies can be attributed in part to the differences in the nuclear matter 
syinmetry energy predicted by non-relativistic and relativistic models. In particular, in Ref. 



Ill ] it has been pointed out that, when the symmetry energy of the RMF models is softened 



to simulate the symmetry energy of Skyrme interactions, a lower value of -ft'nm; consistent 
with the ones used in non-relativistic models, is required to reproduce the energy of the 
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GMR in 208pb. 



Our own results strongly disagree with most of the conclusions of Refs. j9|, llO|, [11 1 . In addi- 
tion to the results obtained in the RMF+RPA framework based on non-linear Lagrangians 
w.h scala. se.f-,Me..ac.o,. B fl Q. in Ref. Q we have cabled ou. ea.cu.a.o„. of 

the isoscalar monopole, isovector dipole and isoscalar quadrupole response of ^°'^Pb, in the 
fully self-consistent relativistic RPA framework based on effective interactions with a phe- 
nomenological density dependence adjusted to nuclear matter and ground-state properties 
of spherical nuclei. The analysis of the isoscalar monopole response with density dependent 
coupling constants has shown that only interactions with the nuclear matter compression 
modulus in the range -ft'nm ~ 260 — 270 MeV, reproduce the experimental excitation energy 
of the GMR in ^''^Pb. In addition, the comparison with the experimental excitation energy 
of the isovector dipole resonance has constrained the volume asymmetry to the interval 34 
MeV < 04 < 36 MeV. 

Even though we cannot explain the difference between the nuclear matter compressibility 
predicted by non-relativistic and relativistic mean-field plus RPA calculations, at least we 
can clarify the apparent inconsistency of the results obtained in the relativistic framework. 
The purpose of this work is twofold. By using the RMF+RPA based on effective interac- 
tions with density- dependent meson-nucleon couplings, we will show: (1) contrary to the 



claims made in Ref. 



our knowledge of the symmetry energy at low density is not that 



poor, and the volume asymmetry cannot be lowered to the range of values < 30 MeV, 
for which relativistic models with < 230 MeV would reproduce the excitation energy 
of the GMR in ^ospb and, (2) relativistic mean-field effective interactions adjusted to repro- 
duce ground-state properties of spherical nuclei (binding energies, charge radii, differences 
between neutron and proton radii), cannot reproduce the excitation energies of GMR if 
Knui < 250 MeV. Therefore, we will reinforce our result that -ft'nm = 250 MeV is the lower 
bound for the nuclear matter compression modulus in nuclear structure models based on 
the relativistic mean-field approximation. 
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II. EFFECTIVE INTERACTIONS WITH DENSITY-DEPENDENT MESON- 
NUCLEON COUPLINGS AND THE RELATIVISTIC RPA 



The relativistic random phase approximation (RRPA) will be used to calculate the 
isoscalar monopole and isovector dipole strength distributions in ^°*Pb. The RRPA rep- 
resents the small amplitude limit of the time-dependent relativistic mean-field theory. A 
self-consistent calculations ensures that the same correlations which define the ground-state 
properties, also determine the behavior of small deviations from the equilibrium. The same 
effective Lagrangian generates the Dirac-Hartree single-particle spectrum and the residual 
particle-hole interaction. In Ref. it has been shown that an RRPA calculation, consistent 
with the mean-field model in the no — sea approximation, necessitates configuration spaces 
that include both particle-hole pairs and pairs formed from occupied states and negative- 
energy states. The contributions from configurations built from occupied positive-energy 
states and negative-energy states are essential for current conservation and the decoupling 
of the spurious state. In addition, configurations which include negative-energy states give 
an important contribution to the collectivity of excited states. In two recent studies P. llJ] 
we have shown that a fully consistent inclusion of the Dirac sea of negative energy states in 
the RRPA is crucial for a quantitative comparison with the experimental excitation energies 
of isoscalar giant resonances. 

The second requisite for a successful application of the RRPA in the description of dy- 
namical properties of nuclei is the use of effective Lagrangians with non-linear meson self- 
interactions, or Lagrangians characterized by density- dependent meson-nucleon vertex func- 
tions. Even though several RRPA implementations have been available for almost twenty 
years, techniques which enable the inclusion of non-linear meson interaction terms in the 
RRPA have been developed only recently 3, 1^ Q| • In Ref. the RRPA matrix equa- 
tions have been derived for an effective Lagrangian with density-dependent meson-nucleon 
couplings. 

Already in Ref. P] we used Lagrangians with non-linear meson self-interactions in time- 
dependent RMF calculations of monopole oscillations of spherical nuclei. The energies of 
the GMR were determined from the Fourier power spectra of the time-dependent isoscalar 
monopole moments < > (t). It has been shown that the GMR in heavy nuclei, as 
well as the empirical excitation energy curve ^ 80 A~^^^ MeV, are best reproduced 
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by an effective force with K„t^ ^ 250 — 270 MeV. This result has been confirmed by the 
RRPA calculation of Ref. Isf. In particular, the best results have been obtained with 

n 

the NL3 effective interaction [16[ (-ft'nm = 272 MeV). Many calculations of ground-state 
properties and excited states, performed by different groups, have shown that NL3 is the 
best non-linear relativistic effective interaction so far, both for nuclei at and away from 
the line of /^-stability. On the other hand, by using the same type of effective Lagrangians 
with non-linear isoscalar-scalar interaction terms, in the recent RRPA analysis of nuclear 
compressional modes of Refs. 0, Q it has been suggested that models of nuclear structure 
having -ft'nm well above ~ 200 MeV are likely to be in conflict with experiment. 

Even though models with isoscalar-scalar meson self-interactions have been used in most 
applications of RMF to nuclear structure in the last fifteen years, they present well known 
limitations. In the isovector channel, for instance, these interactions are characterized by 
large values of the symmetry energy at saturation - 37.9 MeV for NL3 - as compared with 
the empirical value 04 = 30 ± 4 MeV. This is because the isovector channel of these effective 
forces is parameterized by a single constant: the p-meson nucleon coupling Qp. With a 
single parameter in the isovector channel, it is not possible to reproduce simultaneously the 
empirical value of 04 and the masses oi N ^ Z nuclei. On the other hand, extensions of the 
RMF model that include additional interaction terms in the isoscalar and/or the isovector 
channels were not very successful, at least in nuclear structure applications. The reason is 
that the empirical data set of bulk and single-particle properties of finite nuclei can only 
constrain six or seven parameters in the general expansion of an effective Lagrangian (l7l |. 
Adding more interaction terms does not improve the description of finite nuclei. Rather, 
their coupling parameters and even their forms cannot be accurately determined. 

Models based on an effective hadron field theory with medium dependent meson-nucleon 
vertices jl^ ll^ present a very successful alternative to the use of nonlinear self- 
interactions. Such an approach retains the basic structure of the relativistic mean-field 
framework, but could be more directly related to the underlying microscopic description 
of nuclear interactions. In Ref. ^ we have extended the relativistic Hartree-Bogoliubov 



fRHB) model 



22l | to include density dependent meson-nucleon couplings. The effective 



Lagrangian is characterized by a phenomenological density dependence of the cr, uj and p 
meson-nucleon vertex functions, adjusted to properties of nuclear matter and finite nuclei. 
It has been shown that, in comparison with standard RMF effective interactions with non- 
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linear meson-exchange terms, the density-dependent meson-nucleon couphngs significantly 
improve the description of symmetric and asymmetric nuclear matter, and of ground-state 
properties oi N ^ Z nuclei. 

The RRPA with density- dependent meson-nucleon couplings has been derived in Ref . . 
Just as in the static case the single-nucleon Dirac equation includes the additional rearrange- 
ment self-energies that result from the variation of the vertex functionals with respect to 
the nucleon field operators, the explicit density dependence of the meson-nucleon couplings 
introduces rearrangement terms in the residual two-body interaction. Their contribution is 
essential for a quantitative description of excited states. By constructing families of inter- 
actions with some given characteristic (compressibility, symmetry energy, effective mass), it 
has been shown how the comparison of the RRPA results on multipole giant resonances with 
experimental data, can be used to constrain the parameters that characterize the isoscalar 
and isovector channel of the density- dependent effective interactions. In particular, the GMR 
in ^°^Pb requires the compression modulus to be in the range -ft"nm ~ 260 — 270 MeV, and 
the isovector GDR in ^osp]^ jg Q^ly reproduced with the volume asymmetry in the interval 
34 MeV < 04 < 36 MeV. We have not, however, analyzed the influence of the symmetry 
energy on the range of allowed values of -R'nm- We have also not tried to correlate the RRPA 
results for the isovector GDR with data on neutron radii, although it has been shown in 
Ref. ... RHB ca.cu.a«o„. w.h *e den.ty-depe„de. effective interact,™ DD-ME: 
(a4 = 33.1 MeV) reproduce the available data on differences between neutron and proton 
radii for ^osp]^ ^^^^ several Sn isotopes. 



III. NUCLEAR MATTER COMPRESSIBILITY AND SYMMETRY ENERGY 

By performing fully consistent RMF plus RRPA calculations of nuclear ground-state 
properties and excitation energies of giant resonances, in this section we will try to corre- 
late the nuclear matter symmetry energy and the nuclear matter compression modulus of 
relativistic mean-field effective interactions. 



In Ref. 



Piekarewicz has used the standard RMF effective interactions with isoscalar- 
scalar meson self-interactions, and with compression moduli in the range i^'nm ~ 200 — 300 
MeV, to compute the distribution of isoscalar monopole strength in ^°'^Pb. The main result 
of his analysis is that, when the symmetry energy is artificially softened, in an attempt to 
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simulate the symmetry energy of Skyrme interactions, a lower value for the compression 
modulus is obtained, consistent with the predictions of non-relativistic Hartree-Fock plus 
RPA calculations. 

As we have already emphasized in the previous section, the RMF models with isoscalar- 
scalar meson self-interactions are characterized by large values of the symmetry energy at 
saturation volume asymmetry. When adjusting such an effective interaction to properties 
of nuclear matter and bulk ground-state properties of finite nuclei (binding energy, charge 
radius) it would be clearly desirable to bring down the value of to its empirical value = 
30 ±4 MeV. This, however, is not feasible. The simple reason is that, if is brought below ~ 
36 — 37 MeV by simply reducing the single coupling constant Qp in the isovector channel, then 
it is no longer possible to reproduce the relative positions of the neutron and proton Fermi 
levels in finite nuclei, i.e. the calculated masses of N Z nuclei display large deviations from 
the experimental values. The binding energies are only reproduced if a density dependence 
is included in the p-meson coupling, or a non-linear p-meson self-interaction is included in 
the model. If the interaction, however, is adjusted to a single nucleus, as it was done in 
Ref. for ^°^Pb, then even the standard RMF interactions contain enough parameters to 
obtain almost any combination of symmetry energy and nuclear matter compressibility. 

In the present analysis we have used effective interactions with density-dependent meson- 
nucleon vertex functions. For the density dep endence of the meson-nucleon couplings we 
adopt the functional used in Refs. |20, |2l| . The coupling of the cr-meson and tu-meson 



to the nucleon field reads 

9iip) = 9iips^.t)fi{x) for i = a,uj , (2) 

where 

1 + bjjx + djf 

is a function of x = p/psat? and psat denotes the baryon density at saturation in symmetric 
nuclear matter. The eight real parameters in Q are not independent. The five constraints 
= 1, /"(I) = /"(I), and /j"(0) = 0, reduce the number of independent parameters to 
three. Three additional parameters in the isoscalar channel are: ga{psat), QujiPsat), and rria 
- the mass of the phenomenological sigma-meson. For the p-meson coupling the functional 
form of the density dependence is suggested by Dirac-Brueckner calculations of asymmetric 
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nuclear matter l23| 

9p{p) = ^p(Psat) exp [-ap{x - 1)] . (4) 

The isovector channel is parameterized by gp{psat) and ap. Usually the free values are used 
for the masses of the u and p mesons: m^^ = 783 MeV and nip = 763 MeV. In principle one 
could also consider the density dependence of the meson masses. However, since the effective 
meson-nucleon coupling in nuclear matter is determined by the ratio g/m, the choice of a 
phenomenological density dependence of the couplings makes an explicit density dependence 
of the masses redundant. 

Obviously, the framework of density-dependent interactions is more general than the 
standard RMF models, and it encloses models with non-linear meson self-interactions in 
the isoscalar- scalar, isoscalar-vector, and isovector-vector channels. The eight independent 
parameters, seven coupling parameters and the mass of the a-meson, are adjusted to re- 
produce the properties of symmetric and asymmetric nuclear matter, binding energies and 
charge radii of spherical nuclei. For the open-shell nuclei pairing correlations are treated in 
the BCS approximation with empirical pairing gaps (five-point formula). 

In order to investigate possible correlations between the nuclear matter symmetry energy 
and the compression modulus, we have constructed three families of interactions, with iCnm = 
230, 250, and 270 MeV, respectively. For each value of K^^a we have adjusted five interactions 
with 04 = 30, 32, 34, 36 and 38 MeV, respectively. These interactions have been fitted to 
properties of nuclear matter (the binding energy E/A = 16 MeV (5%), the saturation 
density psat = 0.153 fm^'^ (5%), the compression modulus -ft'nm (0.1%), and the volume 
asymmetry 04 (0.1%)), and to the binding energies (0.1%) and charge radii (0.2%) of ten 
spherical nuclei: ^^Q, ^OQa, ^^Zr, Ii2,ii6,i24,i32gj^^ 204,208,2i4p^^ ^j^^ differences 

between neutron and proton radii (10%) for the nuclei ^^^Sn, ^^^Sn and ^°^Pb. The values in 
parentheses correspond to the error bars used in the fitting procedure. Note that, in order to 
fix an interaction with particular values of -ft'nm and 04, a very small error bar of only 0.1% 
is used for these two quantities. For each family of interactions with a given i^nm and for 
each 04, in Fig. [T]we plot the differences, expressed as a percentage, between the calculated 



and experimental binding energies 2^. The corresponding deviations of charge radii are 
shown in Fig. El The results are rather good. Most deviations of the binding ener gies are 

< 0.2%, and the differences between the calculated and experimental charge radii |25| are 

< 0.3%, except for ^°Ca. It should be emphasized that, by using the standard RMF models 



with only isoscalar-scalar meson self-interactions, it is simply not possible to construct a set 
of interactions with this span of values of and 04 and the same quality of deviations of 
masses and charge radii from experimental values. 

For the three sets of interactions with -ft^nm = 230, 250, and 270 MeV, in Fig. |31 we 
display the corresponding nuclear matter symmetry energy curves for each choice of the 
volume asymmetry 04. The symmetry energy can be parameterized 



S2{p) = a4 + ^ (P - Psat) + (P - Psat)^ + " " " (5) 

Psat -LoPsat 



The parameter po defines the linear density dependence of the symmetry energy, and AKq 
is the correction to the incompressibility. 

In finite nuclei, among other quantities, the symmetry energy directly determines the 
differences between the neutron and the proton radii. In a recent study of neutron radii 
in non-relativistic and covariant mean-field models 
neutron skin and the symmetry energy has been analyzed. In particular, the analysis has 
shown that there is a very strong linear correlation between the neutron skin thickness in 
^°^Pb and the individual parameters that determine the symmetry energy 5*2 (p): CI4, Po 
and AKn. The empirical value of r„ — Tp in 208p|^ (0.20 ± 0.04 fm from proton scattering 



the linear correlation between the 



0, 



data j22[, and 0.19 ± 0.09 fm from the alpha scattering excitation of the isovector giant 
dipole resonance places the following constraints on the values of the parameters of 
the symmetry energy: 04 30 - 34 MeV, 2 Mev/fm^ < Po < 4 Mev/fm^, and -200 MeV 
< AA^o < -50 MeV. 

For the three sets of interactions with -ft^nm = 230, 250, and 270 MeV, in the two lower 
panels of Fig. lUwe plot the coefficients po and AfCo as functions of the volume asymmetry 
04. As we have shown in Ref. in order to reproduce the bulk properties of spherical 
nuclei, larger values of 04 necessitate an increase of po- It is important to note that only 
in the interval 32 MeV < 04 < 36 MeV, both po and /S.Kq are found within the bounds 
determined by the value of r„ — Tp in ^''^Pb. The increase of po with 04 implies a transition 
from a parabolic to an almost linear density dependence of S2 in the density region p < 0.2 
fm~^ (see Fig. This means, in particular, that the increase of the asymmetry energy 
at saturation point will produce an effective decrease of 5*2 below p ^ fm~^. In Refs. 
Q,!^! it has been shown that, as a result of the increase of po with 04, the excitation energy 
of the isovector GDR decreases with increasing S'2(psat) = 0,^, because this increase implies 
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a decrease of 5*2 at low densities characteristic for surface modes. This effect is illustrated in 
the upper panel of Fig. |1[ where we plot the calculated excitation energies of the isovector 
GDR in ^°'^Pb as functions of 04, for each set of interactions with -ft'nm = 230, 250, and 270 
MeV. The calculated centroids of the Lorentzian folded strength distributions are shown in 
comparison with the experimental value of 13.3 ±0.1 [30j. The RRPA excitation energies of 
the isovector GDR decrease linearly with 04, and the experimental value favors, for all three 
families of interactions, the interval 34 MeV < 04 < 36 MeV for the volume asymmetry. 

The results of fully consistent RRPA calculations of the isoscalar monopole response in 
^°^Pb are shown in the upper panel of Fig. El where we plot the excitation energies of the 
GMR for the three families of interactions with = 230, 250, and 270 MeV, respectively, 
as functions of the volume asymmetry 04. The shaded area denotes the experimental value: 
E = 14.17±0.28 MeV ^. For each interaction, in the lower panel we plot the corresponding 
result for the difference of the neutron and proton radius in ^''^Pb, in comparison with, at 
present, the best experimental value: 0.20 ± 0.04 fm from proton scattering data [2^. The 
calculated radii practically do not depend on the compressibility but, of course, display a 
strong linear dependence on the volume asymmetry 04. The comparison with the experi- 
mental estimate limits the possible values of the symmetry energy at saturation to 32 MeV 
< 04 < 36. This interval for 04 is somewhat wider than the range deduced from the isovec- 
tor GDR but, nevertheless, both quantities exclude values 04 < 30 MeV and 04 > 38 MeV. 
Coming back to the GMR (upper panel of Fig. EJ, we notice that only the set of interactions 
with ii'nm = 270 MeV reproduces the experimental excitation energy of the GMR for all 
values of the volume asymmetry 04. With i^nm = 250 MeV, only for the two lowest values 
of 04 the RRPA results for the GMR excitation energy are found within the bounds of the 
experimental value, -ft'nm = 250 MeV is obviously the lower limit, if not even too low, for 
the nuclear matter compression modulus of the relativistic mean-field effective interactions. 
This result is completely in accordance with our previous results obtained with relativistic 

nn 

effective forces with non-linear meson self-interactions p, and with density-dependent 
interactions [12]. The calculation absolutely excludes the set of interactions with -ft'nm = 
230 MeV, for any value of 04. The calculated energies are simply too low compared with 
the experimental value. 

Finally, in Fig. IHlwe plot, for all three families of interactions, the calculated differences 
between neutron and proton radii of Sn isotopes, as functions of the mass number, in com- 
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parison with available experimental data j3l|. Similar to the result obtained for ^°®Pb, the 
calculated values practically do not depend on the nuclear matter compressibility. While all 
the interactions reproduce the isotopic trend of the experimental data, and we also notice 
that the error bars are rather large, nevertheless the comparison excludes values 04 < 30 
MeV and 04 > 38 MeV. 



IV. SUMMARY AND CONCLUSIONS 



In this work we have clarified the apparent inconsistency of the values of the nuclear mat- 
ter compression modulus deduced from relativistic RPA calculations of the giant monopole 
resonance (GMR) in ^°^Pb. By using standard RMF effective Lagrangians with non-linear 
isoscalar scalar meson self-interactions, in a number of studies (time-dependent RMF, rela- 
tivistic RPA) we have shown that the experimental data on GMR in heavy nuclei, as well as 
the empirical excitation energy curve E^^SO A-^/^ MeV, are best reproduced by an effec- 
tive force with K^^ ^ 250 — 270 MeV . B ■ The best results have been obtained with 
the well known effective interaction NL3 [l6| (i^nm = 272 MeV). On the other hand, by using 
the same type of effective Lagrangians in the fully consistent relativistic RPA calculations 
of the isoscalar monopole response, Piekarewicz has concluded that RMF models with 
well above ~ 200 MeV will be in conflict with experiment 0, , and that the discrepancy 
between the nuclear matter compressibility determined from non-relativistic and relativis- 
tic mean-field plus RPA calculations of GMR excitation energies can be attributed in part 
to the differences in the nuclear matter symmetry energy predicted by non-relativistic and 
relativistic models 

It is well known, however, that the standard RMF Lagrangians, with meson self- 
interactions only in the isoscalar channel, are characterized by large values of the symmetry 
energy at saturation (volume asymmetry) 04. In fact, if the effective interaction in the 
isovector channel is parameterized by the single p-meson nucleon coupling constant, it is not 
possible to simultaneously reduce the value of below ^ 36 — 37 MeV, and still reproduce 
the experimental binding energies of N Z nuclei. In this framework it is simply not possi- 
ble to construct effective interactions with the volume asymmetry in the range of empirical 
values 04 = 30 ± 4 MeV. 

In the present analysis we have used effective interactions with density-dependent meson- 
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nucleoli couplings. RMF models based on an effective hadron field theory with medium 

ri 

dependent meson-nucleon vertices [iSj], provide a much better description of syrnmetric and 



asymmetric nuclear matter, and of ground-state properties oi N ^ Z nuclei jl9l l20l l2ll]. In 
order to investigate possible correlations between the volume asymmetry and the nuclear 
matter compression modulus, we have constructed three sets of effective interactions with 
-ft'nm = 230, 250, and 270 MeV, and for each value of i^nm we have adjusted five interactions 
with 04 = 30, 32, 34, 36 and 38 MeV, respectively. The interactions have been fitted to 
reproduce the nuclear matter saturation properties, as well as the ground states of ten 
spherical nuclei. 

By employing the fully consistent RMF plus RRPA framework with density-dependent 
effective interactions, we have computed the isoscalar monopole and the isovector dipole 
response of ^°^Pb, as well as the differences between the neutron and proton radii for ^osp]-, 
and several Sn isotopes. The comparison of the calculated excitation energies with the 
experimental data on the GMR and isovector GDR in ^''^Pb has shown that: (i) only for 
-ft^nm = 270 MeV the RRPA calculation reproduces the experimental excitation energy of the 
GMR for all values of the volume asymmetry 04, (ii) K-^ra = 250 MeV represents the lower 
limit for the nuclear matter compression modulus of the relativistic mean-field effective 
interactions, (iii) the isovector GDR constrains the volume asymmetry to the interval 34 
MeV < 04 < 36 MeV. In comparison with the available experimental data, the calculated 
differences between neutron and proton radii indicate that the volume asymmetry should 
be in the range 32 MeV < 04 < 36 MeV, and reinforce our conclusion that cannot 
be lowered to a range of values for which relativistic models with -ft'nm < 230 MeV would 
reproduce the excitation energy of the GMR in ^''^Pb. The disagreement between the nuclear 
matter compression moduli predicted by non-relativistic and relativistic mean-field plus RPA 
calculations, cannot be explained by the differences in the volume asymmetry of the non- 
relativistic and relativistic mean-field models. 

The present analysis has confirmed our earlier results that the nuclear matter compression 
modulus of structure models based on the relativistic mean-field approximation should be 
restricted to the interval ^ 250 — 270 MeV, or even slightly higher. 

In addition, we have also shown that, for the relativistic mean-field models, the isovector 
GDR and the available data on differences between neutron and proton radii limit the range 
of the nuclear matter symmetry energy at saturation to 32 MeV < 04 < 36 MeV. It appears 
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that the GDR favors the high end of this interval, but we stress the fact that in the present 
analysis we have not taken into account the influence of the effective mass on the calculated 
excitation energy of the GDR. This is, however, an effect which really goes beyond the mean- 
field approximation. Rather, more accurate data on neutron radii in heavy nuclei would 
provide very useful information of the isovector channel of the effective RMF interactions. 
On the other hand, as it has been shown in Ref. Q|, there is no correlation between the 
effective mass and the excitation energy of the OMR. The choice of the effective mass, 
therefore, cannot influence the nuclear matter compressibility extracted from the GMR. 
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FIG. 1: The deviations (in percent) of the theoretical binding energies of ten spherical nuclei, 
calculated with the three families of interactions with /Cnm = 230, 250, and 270 MeV, from the 
empirical values 3|- The legend relates the different symbols to the volume asymmetries of the 
corresponding effective interactions. 



FIG. 2: Same as in Fig. ^ but for the charge radii compared with the experimental values 



FIG. 4: The isovector GDR excitation energy of ^'^^Pb (upper left panel), the parameter pQ of the 
linear density dependence of the nuclear matter asymmetry energy (lower left), and the correction 
to the incompressibility AKq (lower right), as functions of the volume asymmetry 04. The shaded 
area denotes the experimental isovector GD resonance energy 13.3 it 0.1 MeV. The three sets of 
symbols correspond to the families of interactions with i^nm = 230, 250, and 270 MeV. 



FIG. 5: The RRPA excitation energies of the GMR in Pb as functions of the volume asymmetry 
04, calculated for the three sets of interactions with -fTnm = 230, 250, and 270 MeV. The theoretical 
centroids are shown in comparison with the experimental excitation energy of the monopole reso- 

n 

nance: E = 14.17 it 0.28 MeV 4]. In the lower panel the corresponding results for the difference 
of the neutron and proton radius in ^''^Pb, are plotted in comparison with the experimental value 
0.20 ±0.04 fm bl- 



FIG. 6: The calculated differences between neutron and proton radii of Sn isotopes, calculated 
with the three sets of interactions with i^nm = 230, 250, and 270 MeV, and with the values of the 
volume asymmetry 04 = 30, 32, 34, 36 and 38 MeV. The theoretical values are compared with 
experimental data from Ref . ^ | . 



FIG. 3: S2{p) coefficient ((SJ of the quadratic term of the energy per particle of asymmetric nuclear 
matter, calculated with the three families of effective interactions with K^a\ = 230, 250, and 270 
MeV. The legend relates the different curves to the volume asymmetries of the corresponding 
effective interactions. 
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